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Abstract 

It is shown that a family of models associated with the kicked Harper 
model is relevant for cyclotron resonance experiments in an antidot array. 
For this purpose a simplified model for electronic motion in a related model 
system in presence of a magnetic field and an AC electric field is developed. 
In the limit of strong magnetic field it reduces to a model similar to the kicked 
Harper model. This model is studied numerically and is found to be extremely 
sensitive to the strength of the electric field. In particular, as the strength of 
the electric field is varied a metal - insulator transition may be found. The 
experimental conditions required for this transition are discussed. 

PACS numbers: 05.45.+b, 03.65.Sq, 73.23.-b 



Transport phenomena of a two dimensional electron gas embedded in a periodic 
potential (antidot array) attracted much attention in theoretical and experimental 
studies [|l], ||. Investigations of magnetotransport |^ and photo-conductivity @] 
for GaAs heterostructures show that the electronic dynamics is essentially nonlinear 
with possible transition to a chaotic regime [|l], §]. Exploration of some chaotic 
systems illuminates a deep connection between the electron transport phenomena 
and quantum and classical chaos [jl], |]. Specifically, the kicked Harper model 
(KHM) is an example of a system which is chaotic in the classical limit. It has been 
introduced theoretically in the field of quantum chaos because of its interesting 
spectral and transport phenomena. It is defined by the Hamiltonian 

= LhCosP + KHCOsqY,5{t-n). (1) 

n 

The model and its variants appear naturally for the kicked harmonic oscillator |^ . 
It is of specific interest since it does not follow the KAM picture. Its dynamics 
corresponds to kicks combined with rotations of the four-fold symmetry. This system 
exhibits chaotic motion in a region that increases with Lh and K^. The model was 



subject to extensive theoretical studies 0,[^]-|]TU|. For some regimes of parameters 
its spectrum of quasi-energies is similar in nature to the energy spectrum of the 
Harper model |]^, ^ 0, |T^. It exhibits classical and quantum diffusion as well 



as localization and anomalous diffusion and even ballistic motion |p-||TO|. The 
motivation for the explorations of the kicked Harper model in the field of quantum 
chaos was so far mainly theoretical, because of the variety of interesting phenomena 
that were found. Since the system can be modeled approximately by the kicked 
harmonic oscillator it can be realized experimentally [0. Nevertheless, it has been 
shown recently that this system and models associated with KHM may be realized 



in experiments on cyclotron resonance []T4|. It was proposed there that the KHM 



can be used to model some aspects of electronic motion in antidot arrays. In this 
paper we show that a family of models associated with KHM is relevant for cyclotron 
resonance experiments 0] in antidot arrays for some conditions. 

The one particle Hamiltonian describing electrons in external fields and in a periodic 
potential is 

H=^Xv--Hr.t))\v{v)^ (2) 

where p = {px,Py) is the two dimensional momentum of an electron with effec- 
tive mass m* and charge e, A(r, t) = {^ExSinutjxB — ^Ey cos 1/1,0) is the vector 
potential in such a gauge that takes into account a constant magnetic field in the z- 
direction B = (0,0,5) and an alternating electric field E = {—Ej-cosut, EySinut) 
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with a frequency z/, while V^(r) = V{x, y) is a two dimensional periodic potential 
modeling the antidot structure. For simplicity the potential 

V{y) = VxCO?, hVj^cos— (3) 

where a, h are periods in the x and y directions, will be considered, although it should 
be much sharper and contain more Fourier components for the antidot lattice 0-0. 
The equations of motion for @) are 

*( * * I ^z? +\ 9^(^^y) . Px e 
= 00 [Py — uj m X -\ — hy cos ut) , x = sm ut, 

V ox m* m*u 

dV Py e 

= - — , ?/ = — - ^ X + — cos i/t, (4) 
oy m* m*v 

where cu* = ^ is the electronic cyclotron frequency. The following change of 
variables y = y + yo sin ut with yo = is useful. The equations of motion for the 
y component in terms of this variable are 

y = ^-uj X, py = — — . (5) 

m* oy 

The fact that ^ = -§- was used. Analogously, we make change of variables px = 
Px + sin ut, leading to 

dV(x, y) 1 uj*6 6 

Px=uj*{py-xuj*m*) — ^ — , x=—{px^ -Eysinvt ExSinvt). (6) 

Ox m* V 



If the amplitudes 5, Ex^ Ey satisfy 

-Ey = Ex, (7) 



UJ" 



V 



the equation of motion for x is particularly simple and reduces to i; = Px/m* . The 
equations of motion result from the effective Hamiltonian 

T-C = ^^\pl + iPy - ^*m*xY] + Vx cos + Vy cos {-^iy + 2/o sin vt)) 

= noiPx,Py,x) + V{x,y,t). (8) 

The tilde is omitted in (^) and in what follows for the sake of simplicity of nota- 
tion. Here Ho = 'Ho{px,Py,x) corresponds to the integrable system describing the 
cyclotron motion, while V = V{x, y, t) is the perturbation leading to classical chaos. 
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As the Hamiltonian 7i is periodic in time with the period T = ^ the following 
quantum mechanical analysis will be carried out in the framework of the Floquet 
theory. The eigenvalue problem of the Floquet operator 

F = -ih^ + H, + V = F^ + V (9) 



is considered. An unperturbed basis of Fq is 



|j,s,M>=^e-'J'^*|s,M>. (10) 
V ^vr 



The wave function > in the coordinate representation is [|T2l ITH]: 



< xMs,u>= ^e^^«/V.(f - y), (11) 

where / = ^hc/eB is the magnetic length, while ipsiz) = [exp(— 2:^/2) /y^y^^^^^-f^sl^) 
is a parabolic cylinder function and Hs{z) is the s-th Hermite polynomial. The ma- 
trix of Fo in this basis is diagonal, 

< f , s' ,u\FQ\j, s,u >= [hu*{s + 2^ ~ ^^j]^s,s' ■ ■ 5{u — u). (12) 

The matrix elements of V can be found with the help of the generating function of 
the Hermite polynomials and take the form 

< j', s', u'|V"|j, s,M >= VxP^giha'^) cos{hu — — |s' — s|) ■ 5jj> ■ 5{u — u') 

+^[P^,s{^0i)J.j..y{K)5{u -u' + 271) + {-iy'-'P;;Xha)Jj,^j{K)5{u -u' - 27r)].(13) 

where Jn{,z) is the Bessel function, and 

P±^(w) = e-™/2^(2vrw)^^Ll''-«'l(±7rtt;) (14) 

for s' > s, while for s' < s one should interchange s' and s. Here L^{z) are the 
associated Laguerre polynomials, while a = b/a and the effective Planck constant 
is h = = The magnetic flux quantum is $o = " and the flux through a 
unit cell of the periodic structure is $ = abB. Because of the periodicity in y the 
allowed transitions change u only by integer multiples of 27r. This is also clear from 



3 



( p!3D . For this reason it is convenient to express u in the form u = 27rn + i) where n 
is an integer while < i} < 27r. During the evolution is constant. 



The eigenstates \\,{} > of the Floquet operator F satisfying F\\,{} >= XIX,!^ > are 
decomposed into the unperturbed basis states 

\K^>=T.(^n,s\j,s,n,i^> . (15) 

Projecting on the state < j,s,n,'d\ one finds that the coefficients C;^ ., satisfy the 
following equation: 

[a + i^M^s = LY.KrCha'') cos [hi27Tn + ^) - - r|]<, 

+ f T.[P^Ahc^)Jf-ji'^yn+l,r + {-ir-'-p-r{ha)J,.,,{K)ClA (16) 

The following change of variables was used here: v/uj* hv and h\= [\ — huj*{s + 
|)]/[/;.ti;*], while L = , /C, and K = Jy- ■ /C. The parameter lilC = JV^Vy/hu* is 



the strength of the coupling of the Landau bands introduced in for a stationary 



problem of Hall conductance in this system. In the case when it is not too large 



(/C < 6) the coupling between the Landau bands is weak |jT2[. In the quantum case 
when ha < 1 and ha~^ < 1 we can neglect this interaction, i.e. the terms with 
s 7^ s' in (p!6D can be omitted. This can be seen from (p!4D since the leading term of 



Ps',siui) behaves as w^^' Within this approximation (1I6|) reduces to: 

[A + uj]c^^ = Lcos[h{27rn + + ^Yl J,-/(«:)[(-l)^-^'c{+i + (17) 

j' 

where L = LP+^{ha''^) / P^.,iha), K = K, and [A + uj]/ p-,.(ha) ^ [A + uj], while 
the index s was suppressed in the notation. The Hamiltonian 

^1 = L cosp + _ft'cos(g + Ksin i/if:), (18) 

reduces to (0) in the basis e^'^'^e"*-'^*, where k = ^ = ^^^^^ and the canonical 
variables are q and p = —i27ih^ satisfying [p, q] = —i2'nh. In the absence of the AC 
electric field k = and (|1^) reduces to the Hamiltonian of the Harper model. For the 
model of this paper the semi-classical limit /i — > is the limit of strong magnetic field 
$ — i> cxD (or 5 — > oo). The validity of equations ( p!7D and (0) requires ^ = /i < 1, 
therefore the classical limit is meaningful for (|18]). This Hamiltonian was studied in 
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1^. It was derived there for a model opposite to the one studied here, namely the 
tight binding model where the external fields can be considered as a perturbation 
on the lattice potential. 

The driving potential K cos{q — Ksinut) (related to the one of (|TH|) by a trivial 
shifting of time) is well known in the literature and it has been discussed in the 
context of the description of dynamical localization in atomic momentum transfer 



1^, |T^. This effect has been observed experimentally following an extension of a 
theoretical proposal [Q. The Hamilton equations for Hi of ([ISD are q = —Lsinp 
and p = Ksimplt) with ip{t) = q — Ksinut. For k ^ 1 the forcing resulting in 



change of momentum is dominated by the resonant points [0, where = or 

— Lsinp = q = Ku cos ut. (19) 

Expanding ilj{t) around the resonant point and integrating the Hamilton equa- 
tion for p, taking into account the fact that this integral accumulates most of its 
contribution from a narrow region (of width |/tz/^ sin z/t^l"^^^) around the resonance, 
one finds that the momentum transferred at each resonance is 



^Pr = \ 21^-"" , \ Ksin{^f ± ^), (20) 

y hiv^ I sm vtr I 4 
where V'r = '^{tr) and the sign depends on the direction of crossing of the resonance. 

The position of the resonance depends on p. For strong driving so that ku ^ L for 
the resonance condition it is required that cos utr ~ or | sinz/t^l ~ 1- Consequently 
the resonant points are ~ ±f and therefore are approximately equally spaced 
in time, and occur at the times: ut' = — f + 27r/ and ut^ = f + ^tt/ (/ are integers). 
The resulting map is M = M2 ■ Mi with 



Ml 



p + Ki sin(g + Kq) 
q — Li sinpi 



. { P2= Pi + Ki sin(gi - kq) ,^1) 

^ ■ I 52 = gi - ^isinpa 

where Li = while Ki = ^J^^K and = k — f . First we note that the map is 
periodic in k with the period 2tt and therefore it is periodic in the magnitude of the 
electric field. This map is generated by the Hamiltonian: 

r °° TT 
0.2 = L cosp + pKi ■ [ cos(g + kq) ^ diyt H 27rn) 

n=— 00 ^ 

00 ^ 

+ cos(g-Ko) 5{ut- --2Tm)\. (22) 
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In the case when kq = this system corresponds to the KHM. Numerical analysis 
shows that the model is extremely sensitive to the parameter hq in both quantum 
[h ~ 1) and classical {h = 0) regimes. In the quantum regime we study numerically 
the time evolution of the variance 

var{t) =< n\t) >-< n{t) >'= E^'l/nWP (23) 

n 

and energy spreading over the unperturbed level spectrum as a result of the evo- 
lution from the initial level occupation /„(0) in the momentum space. The stan- 
dard technique of the fast Fourier transform is used to evolve the quantum map 
( |22D of the wave function over the period T over which two kicks take place: 
'if{t + T) = f/\E'(t), where U is the evolution operator over the period and \I'(t, q) = 
S^-oo 6^™''/n(^)- The initial distribution used is /n(0) = 6nfl- An essential 
deviation from KHM is found as the parameter k changes. The parameters Ki, Li 
and h are the same as in p and are explicitly specified in the figures. 

It is found here that when the parameter kq is varied, localization-delocalization 
transition takes place, with different realizations of the delocalization as diffusion, 
anomalous diffusion and ballistic motion. Such a transition from diffusive motion 
for = to ballistic one for kq = gn where g = golden mean is shown 

in Fig. 1. A localization-delocalization transition is found also for Ki/Li < 2 { see 
Fig. 2), such that for KHM {kq = 0) there is dynamical localization of quasi-energies, 
while ballistic motion is found for kq = gn. When Ki,Li -C 1 a slow diffusive 
process in momentum (reported first in Ref. 0) is observed for KHM in the semi- 
classical limit. It is suppressed when the parameter kq increases from zero as shown 
in Fig. 3 resulting in transport motion in real space. This transition corresponds 
to a metal-insulator transition in the |/-direction of the sample. Note the relation 
between the components of the electric field that is crucial for the preference 
of the y direction. It is enhanced in the limit when h = 0. Numerical iterations 
of the classical map ( pT]) for small values Ki and Li with initial distribution taken 
in the vicinity of a separatrix are presented in Figs 4 and 5. These calculations 
demonstrate that this transport in the y direction found for finite values of and 
absent when it vanishes (for the same values of Ki and Li) is of pure classical origin. 
Therefore it is expected to be robust against effects of noise. 

A model for transport in a two dimensional electron gas embedded in the periodic 
super-lattice potential in the presence of external fields was studied. The motion cor- 
responds to cyclotron resonance dynamics with possible metal-insulator transitions 
induced by the high frequency driving. It can be realized experimentally by a 2D 
periodic andidot array. The required conditions for experimental realization of this 
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effect can be acliieved [0. Form* ~ O.lme, a, 6 ~ 10^ ^4- lO^nm, B ~ 0.07^0.37T, 
one obtains tliat uj* ~ 10" ^ 5 • 10"s"\ wliile h ~ 0.08 6. For u ~ lO^s"^ and 
Eq ~ ICGSE one obtains tliat k ~ 300 and tlie condition k ^ 1 is fulfilled. The 

amplitudes of the periodic potential can satisfy the condition /C = '^l^j^'^^ - < 6 
for Vx, Vy < 4meV and a ^ b lO^nm. Reasonable values for the parameters of 
( ^11) and (p2D are Ki = Li = 5 that are obtained for Vx = O.lSmeV and Vy = 3meV 
for a = 6 ~ lO^nm, = 300 and u = lO^^s"^. These were used in numerical 
calculations of the present paper. The values of Ki and Li can be reduced and also 
numerical calculations for smaller values were performed. 

Two metal-insulator transitions induced by the variation of Kq are predicted. The 
first one is a quantum effect of strong delocalization in p corresponding to anomalous 
(super) diffusion or acceleration for p, so that < {p — PoY >~ with 1 < fx < 2 
for 7^ and for any ratio between Ki and Li, while it is known that for kq = 0, 
that corresponds to the KHM, super-diffusion takes place only for Ki > Li. It 
means that for 7^ an initial wave packet spreads in the direction of stronger 
modulation of the 2D super-lattice potential. So far, in the one-band structure 
systems both integrable like the Harper model and chaotic like the kicked Harper 
model (kq = 0), wave packets were found to spread in the direction of the weaker 
modulation of the potential amplitudes. Recently such effect has been studied in 
the many-band structure system with strong interaction between the bands leading 
to chaos in the classical limit |]T^. It corresponds to electronic transport in a 2D 



periodic potential with different modulation amplitudes in x- and y- directions that 
are perpendicular to the magnetic field. The second one is a classical effect taking 
place for i^i,Li < 1 in the both classical h = and semi-classical limit /i <^ 1. It 
corresponds to the strong delocalization in the q direction with anomalous diffusion 
or ballistic motion at any ratio for kq 7^ while for kq = anomalous diffusion 
in q takes place only when < 1. This effect can be explained by the topological 
reconstruction of phase space due to a variation of leading to different settling 
of elliptic and hyperbolic points. An analogous phenomenon has been studied for 
the web map |^ and it has been shown that there is a connection between the 



web structure and features of diffusion in phase space [|T0|- For the first case the 
quantum transition due to the non- vanishing kq corresponds to delocalization in the 
p direction and localization in q. The second transition corresponds to delocalization 
in q and localization in p. Some caution is required when the model studied here 
is compared to antidot lattices, where the potential ^(r) in (0) should be replaced 
by a sharper function. In this case the equation ( |T7|) and the effective Hamiltonian 
(II8D have more complicated form due to interaction between x and y degrees of 



freedom. For a sufficiently strong magnetic field, so that ha and ha ^ are much 
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smaller than unity, and /C < 6 the one Landau band approximation is valid 
and the metal-insulator transitions considered here can be observed. (This justifies 
ignoring the terms with s 7^ s' in ([T3|).) The result may be modified since the form 
of the equation corresponding to (|1^) may be somewhat more complicated because 
of the sharpness of the potential. 
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Figure Captions 



Fig. 1 The quantum occupation probability \fn{t)\ vs n and the variance var{t) 
found from integration of (|22|) for Ki = Li = 5 and h = with = (a,b) 
and Hq = ng (c,d). 

Fig. 2 Same as Fig. 1 but for Li = 4 and K2 = 2. 

Fig. 3 Same as Fig. 1 for i^i = Li = | abd h = for = (a,b) and = gn 
(c,d). 

Fig. 4 Classical evolution of the map (^) for Ki = Li = ^ and = 0. Eleven 
trajectories are presented in figure (a) and folded to the 27r x 27r torus in figure 
(b). The variances are presented in (c) and (d). 

Fig. 5 Same as Fig. 4 but for kq = 0.8 ■ gn. 
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